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NUMERICAL  ME!  HODS  IN  COMPUTATION  OF  RADAR 
DEI  ECTION  PERFORMANCE 
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INTRODUCTION 

The  fornuilas  to  compute  radar  detection  performance  can  be  found  in  many  literatures, 
such  as  [I  J-[9],  Ref.  ( I]  uses  those  formulas  (o  plot  a  large  amount  of  figures  which  show  the 
probability  of  detection  vs.  signal-to-noise  ratio  for  various  target  models,  pulse  numbers,  and 
cell  numbers.  These  figures  are  very  important  in  design  and  development  of  radar  systems, 
as  well  as  in  radar  signal  processing.  The  target  models  studied  in  Ref.  1 1 J  include 

1)  Single  pulse,  con'^lant  amplitude  scattcrers, 

2)  Single  pulse,  Rayleigh  scattcrers, 

.3)  Single  pulse,  dominant  plus  Rayleigh  scattcrers, 

4)  Multi-pulses,  constant  amplitude  scattcrers, 

.5)  Multi-pulses,  slow  fluctuating  targets,  Rayleigh  scattcrers, 

6)  Multi-pulses,  fast  fluctuating  targets,  Rayleigh  scatlercrs, 

7)  Multi-pulses,  slow  fluctuating  targets,  dominant  plus  Rayleigh  scattcrers, 

8)  Multi-pulses,  fast  fluctuating  targets,  dominant  plus  Rayleigh  scattcrers. 

A  serious  difficulty  existing  in  radar  detection  performance  calculation  is  that  many 
formulas  involved  have  very  bad  numerical  behaviors,  and  thus  lead  to  significant  errors  of  the 
results.  Typical  numerical  error  sources  arc  truncation  of  numbers,  subtraction  of  nearly  equal 
numbers,  truncation  of  scries,  numerical  integration,  overflow  underflow,  and  open-range 
integration.  The  overflows  and  underflows  are  mainly  caused  by  two  involved  .special 
functions — incomplete  Tlironto  functions  and  modified  Bessel  functions. 

Most  computers  can  handle  overflows  and  underflows  automatically.  The  computers 
assign  a  maximum  representable  value  to  the  variable  which  generates  overflow,  and  assign 
zero  to  the  valuable  which  generates  underflow.  In  most  ca.ses,  overflow  will  lead  to  significant 
error  of  the  results.  And  in  some  cases,  underflow  will  lead  to  significant  errors.  We  can 
identify  the  underflow  ca.ses  which  will  lead  to  error  and  the  cases  which  will  not.  T)  ensure 
(he  accuracy  of  the  results,  we  must  eliminate  all  the  overflow  cases  and  the  underflow  cases 
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which  will  lead  to  error. 

In  this  report,  wc  will  analyze  the  causes  of  the  numerical  dirficullics,  and  develop  the 
methods  to  overcome  these  difficulties.  The  methods  are  used  to  compute  (he  probability  of 
detection  for  various  target  models,  pulse  numbers,  and  cell  numbers.  Many  of  the  figures  in 
this  report  are  plotted  to  show  the  probability  of  detection  against  the  signal-to-noise  ratio. 
The  figures  in  this  report  are  compared  with  those  from  (I).  Many  of  the  figures  show 
significant  difference  from  those  of  [I  j.  One  of  the  pur}>oses  of  the  work  in  this  part  of  the 
report  is  to  use  the  numerical  methods  in  Part  II  of  this  report,  which  investigates  the  effects 
of  energy  distribution  on  the  radar  detection  performance.  The  methods  developed  here  can 
also  be  applied  to  the  problems  other  than  radar  detection. 


niE  INVOLVED  FORMUIA 

The  formula  involved  in  compulation  of  radar  detection  performance  can  be  found  In 
various  literatures.  The  formulas  listed  below  arc  obtained  from  Ref.  ( 1 ).  It  is  assumed  that  the 
energy  is  uniformly  distributed  over  the  resolution  cells. 

The  relatioaship  between  false  alarm  probability,  Pfa,  and  threshold,  Yi„  is  given  by 

for  single  pulse  cases  (II  j,) 

for  multi-pulse  cases  (lib) 


where  M  is  the  total  number  of  resolution  cells;  N  is  the  number  of  pulses.  l(u,  s)  is  the 
incomplete  Gamma  function  expressed  by 
,  u  vTTs” 


I(u,  s)  =  f 
'  0 


e  *  t® 


s! 


dl 


(1.2) 


The  probability  of  detection,  P^^,  is  given  by 

M 

clni^ 


Pcl= 


OJ) 


where  is  Ihe  probabilily  of  a  single  cell  lo  accross  <hc  Ihrcshold,  and  has  differcnl 
expressions  for  different  large!  models,  as  shown  below. 

I)  For  single-pulse  conslani  ampliliidc  scallcrcrs,  Pj,„  is  given  by 

Pcim  =  I  -  T  I  1,0,  n^e/cMN.)  J  (••4) 

I  II 


where  E  is  (he  (olal  signal  energy.  No  is  Ihe  lolal  noise,  and  TH(m,n,rI  is  (he  inconiplclc 
Toronto  function  which  has  the  following  expression; 


T  (m,n,r) 
0 


(1.3) 


and  B„(z)  is  the  nth  order  modified  Bessel  function  expressed  by  cither 

OO 

I  n  +  2k 


B(7.)=  I  _ 

"  k  =  o  k!(k  +  n)! 


(7/25 


(1.6a) 


or 


_  (7/2)" 

■  r(n+  l/2)  r(l/2)  ''-1 


r '  2  11-1/2  -7.1  . 

J  (1-1  )  c  dt 


(1.6b) 


where  r(n)  is  the  Gamma  function. 

2)  For  single  pulse  Rayleigh  scattcrers,  Pj^,  is  given  by 

Y 

Pdm  =  exp(-  ' 


(  1  +  e/(2MN„)) 


(1.7) 


y)  For  single-pulse  dominant  plus  Rayleigh  scattcrers,  Pj,,,  is  given  by 


dm  — 


I  +  4MN(/E 


r( 


4MN, 


1  + 


Y 


I  + 


^(4MNo) 


) 


exp^ 


I  + 


^(4MNo) 


) 


(l.«) 


4)  For  multi-pulse  constant  amplitude  scalterers,  Pjn,  is  given  by 
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(1.9) 


^  '^ne  /(MN^)  I 


where  N  is  (he  number  of  (he  pulses. 

.5)  For  muKi-pulse  slow-nuc(ua(ing-(arge(  Rayleigli  sca((erers,  Pd,,,  is  given  by 

N-1 


1 


NE/(2MN, 

1  +NE/(2MN()))*  *(n/n^ 


Y„ 


■|(I  +  2MN^(NE)) 


erH 


(I- 10) 


6)  For  muKi-pulse  fas(-fluc(ua(ing-(argc(  Rayleigh  seaderers,  is  given  by 

Y 


Pdm  =1-1 


(n/N  (l  +  E/(2MN^)) 


(1.11) 


7)  For  muhi-pulse  slow-nuc(ua(ing-(afgc(  dominan(  plus  Rayleigh  sea((eiers,  P{jni  is 
given  by 


00 


Pdm  =  f  P(Y)dY 
Yh 


where 


P(Y)  = 


/  _ 1 _  ^N-2 

7.rNE-/(4MN^)j  ^  r 

n  +NF  /r4MN V(I  + 


-Y 


+ 


(l+NE/(4MN^y 

(N-2)(l+  NE-/(4MN„)) 
(14-NE/(4MN^)2 
Y^“^  exp(-Y) 


NE/(4MN,j))) 
NE/(4MN^)) 


(N-2)!  (I  +  nE  /(4MN^))' 


(1.12a) 


(1.12b) 


and 
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K=  .(- 


(l  +  4MrV(NE)) 


(l.l2c) 


8)  For  muHi-pulse  fast-fluctuating-largcl  dominanJ  plus  Rayleigh  scallcrcrs,  Pjn,  is 
given  by 


^  dm  —  1 


N! 


N 


(1  +  E/(4Mr^))'^  k  =  0 

'^44  +  K(I  +  E/(': 


E/(4MN^)) 


!  (  4MN„  ) 

) 


,  N  +  k-1 


(113) 


THE  NUMERICAI.  DIFFICULTIES 

The  formulas  involved  in  compulalion  of  radar  detection  arc  given  in  the  above  scclion. 
Most  numerical  difficulties  arc  caused  by  the  incomplete  Tbronto  function,  modified  Bessel 
function  and  the  open-range  integral. 

The  incomplete  Tbronto  function  and  the  Bessel  function  will  generate  many  overflows 
and  underflows.  For  the  open-range  integral  in  Eq.  (I.l2a),  the  function  P(Y)  may  decay 
slowly  with  Y.  Therefore,  in  order  to  obtain  accurate  results,  we  mast  take  a  large  upper  limit 
of  the  integral.  However,  when  the  variable  Y  becomes  large,  underflows  and  overflows  may 
occur  in  computing  P(Y).  Also  large  integral  range  will  result  in  large  computational  time. 

In  the  incomplete  Tbronto  function  (sec  Eq.  (1.5)),  the  following  problems  may  appear; 

1.  The  factor  r"  '  may  underflow  when  r  is  large  and  m  is  much  greater  than  n. 

2.  The  factor  exp(-r^  may  underflow  when  r  is  large. 

3.  In  the  case  where  m  is  much  larger  than  n,  the  factor  f"  "  may  underflow  or  overflow 
when  t  is  small  or  large,  respectively. 

4.  TTk  factor  cxp(-t^  may  underflow  when  t  is  large. 
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In  the  modified  Bessel  funclion,  if  ihc  series  form  (Eq.  (I  da))  is  used  and  z  is  no<  small, 
we  must  sum  up  many  terms  in  order  to  obtain  accurate  results.  Then,  the  factors 
k!,  and  (k  +  n)!  may  overflow.  I f  the  integral  form  (Eq.  ( 1 .6b))  is  used  and  n  is  zero,  there  are 
singular  points  at  t  =  - 1  and  t  =  I .  These  singular  points  will  .seriously  degrade  the  accuracy  of 
the  results. 

In  Eq.  (1. 12b),  when  Y  is  large,  the  factor  '  may  overflow,  the  exponential  factors 
may  underflow.  Also  large  value  of  Y  may  lead  to  problems  in  incomplete  Gamma  function 
(Eq.  1.2)).  In  the  incomplete  Gamma  function,  the  factor  e  ‘  may  underflow  and  the  factor 
may  overflow. 

The  other  numerical  error  sources  include  truncation  of  scries,  truncation  of  numbers, 
and  numerical  integration.  Some  scries,  such  as  the  modified  Be.sscl  function,  may  decay 
slowly  with  the  index.  Then  if  the  scries  arc  truncated  too  early,  the  results  arc  not  accurate; 
while  if  the  series  afc  not  truncated  early,  overflows  or  underflows  may  occur. 

THE  NUMERICAI.  METHODS 

It  is  relatively  easy  to  reduce  the  error  caused  by  truncation  of  numbers  and  numerical 
integral.  We  can  reduce  the  errors  cau.sed  by  truncation  of  numbers  by  using  double  precision 
variables  which  arc  generally  sufficient.  And  we  can  reduce  the  errors  caused  by  numerical 
integration  by  dividing  the  integration  intervals  finely.  Also  if  double  precision  variables  arc 
used,  we  find  that  the  errors  caused  by  subtraction  of  numbers  arc  not  serious. 

To  improve  the  behaviors  of  the  modified  Bessel  function,  we  use  the  integral  form  when 
n  is  not  zero  and  use  the  scries  form  when  n  is  zero.  This  method  can  reduce  large  numbers 
of  overflows  and  underflows  and  to  avoid  the  singular  problems. 

The  remaining  difficulties  arc  overflows  and  underflows.  To  handle  these  problems,  we 
classify  the  overflows  into  the  following  patterns: 

1.  S»t  "  when  t  and  S  arc  small  and  n  is  large.  In  this  case  the  factor  t  "  overflows. 
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Generally,  when  a  faclor  overflows,  (here  should  be  anolher  fae(or  which  is  small,  odicrwise 
(he  problem  can  nol  be  solved.  In  (his  padern  we  can  solve  (he  problem  by  using  an  equivalcnl 

form  of  . 

2.  S*r"  when  rand  n  arc  large,  S  is  small.  In  (his  ease,  (he  fac(or  r"  overflows.  We  can  solve 

(he  problem  by  using  an  equivalcnl  form  of  r"  ^  ‘ . 

.1.  S»X(i  when  some  of  (;  overflows  and  S  is  small.  In  (his  case,  we  can  solve  (he  problem 
by  using  an  cquivaleni  form  of  X(S*(i). 

Similarly,  we  can  classify  underflows  in(o  (he  following  paderns; 

1 .  B«cxp(-  ()  when  B  and  ( arc  large.  In  (his  ca.se,  (he  faclor  cxp(-()  underflows.  Generally 
compulcr  will  assign  zero  (o  (he  undcrflowcd  variables.  This  mclhod  will  nol  inlroducc  errors 
if  no  faclor  is  large.  Therefore,  we  only  need  (o  lake  care  (he  cases  wilh  a  large  faclor,  B.  In 
(his  padern,  we  can  solve  (he  problem  by  using  an  cquivaleni  form  of  cxp(-(  +  /nB).  If  (he 
original  form  is  used,  (he  value  of  (he  expression  is  zero,  biK  (he  cquivaleni  form  may  give 
non-zero  values. 

2.  B*!*^  when  ( is  small,  B  and  n  are  large.  In  (his  case,  (he  faclor  I"  underflows.  We  can 

solve  (he  problem  by  using  an  cquivaleni  form  of  l"  ^  ’ . 

.1.  B»r  "  when  r,  n,  and  B  arc  all  large.  In  (his  case,  (he  faclor  r "  underflows,  and  we  can 

solve  (he  problem  by  using  an  cquivaleni  form  of  r""  ^  . 

4.  B»£(j  when  some  of  (j  underflows  and  (he  sum  is  small.  We  can  solve  (he  problem  by 
using  an  cquivaleni  form  of  X(B«li).  If  (he  sum  is  nol  small,  (he  undcrflowcd  (erm  can  be 
ncglccled  wilhoul  iniroducing  errors. 

By  using  (hesc  mclhod,  we  can  significanlly  improve  (he  accuracy  of  (he  rcsulls.  The 
melhods  can  climinalc  all  (he  overflows.  Allhough  some  underflows  slill  exisi,  (hey  will  nol 
inlroducc  errors;  c.g.,  all  (he  faclors  arc  small,  or  a  Icrm  in  a  scries  is  ncgleclablc  compared 
wilh  olher  (erms  in  (he  series. 
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THE  RESlI»;rS 


There  are  39  independent  figures  for  various  target  models  and  pulse  numbers  in  1 1 1.  The 
figures  shows  the  probability  of  detection  vs.  the  signal  to  noise  ratio  for  different  cell 
numbers.  By  using  the  numerical  methods  developed  above,  we  reproduced  all  of  (hose 
figures. 

We  compared  the  figures  with  (hose  in  1 1 1,  and  found  that  about  36%  of  our  reproduced 
figures  arc  significantly  different  from  (ho.se  in  1 1).  The  others  are  essentially  the  same. 

Figures  1.1-1.14  show  (he  comparison  of  the  figures  with  significant  differences.  Figures 
1.  la-1. 14a  are  the  reproduced  ones  by  using  the  numerical  methods  in  this  report.  Figures 
1.  lb-1. 14b  are  the  corresponding  figures  from  ( 1 ).  Figures  1.15- 1  ..39  show  the  figures  without 
significant  differences  (with  good  numerical  behaviors). 


SUMMARY 

The  formula  involved  in  calculation  of  the  radar  detection  performance  have  bad 
numerical  behaviors,  and  thus  lead  to  errors  of  the  results.  Our  numcriciil  methods  developed 
can  overcome  (hose  difficulties,  and  give  much  more  accurate  results. 

Incomplete  Toronto  function,  modified  Bessel  function  and  open-range  integration  are 
main  sources  of  numerical  problems. 

The  target  models  with  good  numerical  behaviors  includes  single-pulse  Rayleigh 
scatterers  and  single-pulse  dominant  plus  Rayleigh  scaltcrcrs.  The  target  models  with  bad 
numerical  behaviors  include  single-pulse  constant  amplitude  scatterers  and  multi-pulse 
constant  amplitude  scatterers.  The  target  models  with  some  numerical  problems,  especially 
when  N  and  Yb  are  large,  include  multi-pulse  Rayleigh  scaltcrcrs  (slow  fluctuating  or  fast 
fluctuating  target.s)  and  multi-pulse  dominant  plus  Rayleigh  scatterers  (slow  fluctuating  or 
fast  fluctuation  targets). 
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Figure  1.2a  The  reproduced  figure  for  single-pulse  conslanl  amplitude  scalterers 
when  Pfa=  l.E-5 


Figure  1.2b  The  figure  from  Ref.  { 1  j  for  single-pulse  constant  amplitude  scalterers 
when  Pfa  =  l.E-5 


Figure  1.3a  The  reproduced  figure  for  single-pulse  conslanl  amplitude  scallerers 
when  Pfa=  l.E-8 
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Figure  1.3b  The  figure  from  Ref.  f  I )  for  single-pulse  conslanl  amplitude  scallerers 
when  Pfa=  l.E-8 
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Figure  1.4a  The  reprcxluccx!  figure  for  4-piiIse  conslanl  am|)litucle  scatterers 
when  Pfa=  l.E-2 
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Figure  1.5a  The  reprocliiceci  ligure  (or  4-pulse  constant  amplitude  scatterers 
when  Pra=  l.E-5 


rnM»nt«r  •f  MtKtloai  n.  S««Mt  t>  MoIm  K*t1« 


St9n«t  tP  NoHc  lUtto.  I/K.  ei 


Figure  1.5b  The  figure  from  Ref.  1 1 1  for  4-pul.se  constant  amplitude  scatterers 
when  Pfa=  l.E-5 


14 


Figure  1.6a  The  reprotluced  figure  for  4-pulse  constant  amplitude  scatterers 
when  Pfa  =  1  .E-8 
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Figure  1.6b  The  figure  from  Ref.  1 1 1  for  4-pulse  constant  amplitude  scatterers 
when  Pfa=  I.E-<S 
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Figure  1.7a  The  reproduced  figure  for  16-pulse  constant  amplitude  scatterers 
when  Pfa  =  1  .E-2 
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Figure  1.7b  The  figure  from  Ref.  (1 1  for  16-puIse  constant  amplitude  .scatterers 
when  Pfa  =  1  .E-2 
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Figure  l.8a  The  rcprotlucctl  figure  for  16-pul.sc  conslanl  amplilucic  scatlerers 
when  Pfj,  =  1  .E-5 
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Figure  1.8b  The  figure  from  Ref.  ( 1  j  for  16-pulsc  conslanl  amplitude  scallcrers 
when  Pfj,=  l.E-5 
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Figure  1.9a  The  reproduced  figure  for  16-pulse  constant  amplitude  scatterers 
when  Pfa  =  1  .E-8 


Figure  1,9b  The  figure  from  Ref.  [1 1  for  16-pulse  constant  amplitude  scatterers 
when  Pfa=  l.E-8 
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Figure  1.10a  The  reproclucetl  figure  for  I6-pul.se  slow-nucliiating  Rayleigh 
scatlercrs  when  IVa  =  l  .E-8 
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Figure  1.10b  The  figure  from  Ref.  1 1 1  for  16-pulsc  slow-nuclualing  Rayleigh 
scaltcrcrs  when  Pfa  =  I  .E-8 
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Figure  1.11a  The  reproduced  figure  for  16-pulse  fasl-lluclualing  Rayleigh 
scatterers  when  Pfj,  =  l.E-8 
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Figure  1.11b  The  figure  from  Ref.  [1 1  for  l6-pul.se  fasl-lluctuating  Rayleigh 
scatterers  when  Pfj,  =  l.E-S 
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Figure  L12a  The  reproduced  figure  lor  16-pulse  slow-Huctuating  dominant  pi 
Rayleigh  scatterers  when  Pfj,  =  l.E-2 
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Figure  l.l2bThe  figure  from  Ref.  f  1]  for  16-pulscslow-nLiclualing  dominant  pi 
Rayleigh  scatterers  when  Piv,  =  l.E-2 


Figure  1.13a  The  reproduced  figure  for  16-piilseslow-llucluating  dominant  plus 
Rayleigh  scatterers  when  Pf^  =  l.E-8 
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Figure  1.13b  The  figure  from  Ref.  (1|  for  l6-pul.se  slow-nucliiating  dominant  plus 
Rayleigh  .scatterers  when  Pfj,  =  l.E-8 
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Figure  1.14a  The  reproduced  figure  for  16-pulse  fasl-nuclualing  dominant  plus 
Rayleigh  scatterers  when  Pf^  =  l.E-8 
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Figure  1.14b  The  figure  from  Ref.  [1 )  for  16-pulsc  fast-Huclualing  dominant  plus 
Rayleigh  scatterers  when  Pra=  l.E-8 


Figure  1.15  The  reproduced  figure  for  .single-pulse  Rayleigh  scatlcrcrs 
when  Pf{,=  l.E-2 


Figure  1.16  The  reproduced  figure  for  single-pulse  Rayleigh  scallercrs 
when  Pfj,  =  1  .E-5 
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Figure  1.18  The  reprocluceci  figure  for  single-pulse  th^iiiinnnl  plus  Rnyleigh 
scalterers  when  P(;,  =  I  .E-2 
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Figure  1.19  The  rcprotlucccl  figure  for  single-pulse  deminnnl  plus  Rayleigh 
scatterers  when  Pf;,=  l.E-5 
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Figure  1.20  The  repixulucecl  figure  for  single-pulse  dominanl  plus  Rayleigh 
scatterers  when  Pf;,=  l.E-8 
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Figure  1.21  The  reproduced  figure  for  4-pulse  slow-fluclualing  Rayleigh 
scatterers  when  Pf.,=  I.E-2 
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Figure  1.22  The  reprcHlucedi  figure  for  4-puIse  slow-fluctuating  Rayleigh 
scatterers  when  Pr;,=  I.E-5 
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Figure  1.23  The  repnxlucecl  figure  for  4-pulse  slow-lliicliinling  Rnyleigh 
scallerers  when  Pf.,  =  I  .E-S 
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Figure  1.24  The  reproduced  figure  for  16-pulse  slow-lluctualing  Rayleigh 
scallerers  when  Pfj,==  I.E-2 
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Figure  1.25  The  reproduced  figure  lor  16-pulse  slow-fiiictuating  Rayleigh 
scatterers  when  Pf.,=  I.E-5 
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Figure  1.26  The  reproduced  figure  for  4-pulsc  rasl-Hucluating  Rayleigh 
scallerers  when  Pr-,  =  I.E-2 


Figure  1.27  The  reproduced  figure  for  4-pulse  fasl-lluci-iKUing  Rayleigh 
sea tterers  when  Pfj,=  I.E-5 


Figure  1.28  The  reproduced  figure  for4-puIse  fasl-lluclualing  Rayleigh 
sea  tterers  when  Pn,=  I.E-8 
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Figure  1.29  The  reproduced  figure  for  16-pulsc  fast-nucluating  Rayleigh 
scatterers  when  Pf.,  =  1  .E-2 


Figure  1.30  The  reprotluced  ligure  for  16— pulse  fast— nuclualinc  Rayleitih 
scatterers  when  Pf.,  =  I  .E-5 
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Figure  1.31  The  reproduced  figure  for  4-piilse  slow-llucliuning  dominant  plus 
Rayleigh  scalterers  when  Pfj,  =  l.E-2 
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Figure  1.32  The  reproduced  figure  for  4-pul.se  slow-flucliialing  dominant  plus 
Rayleigh  scatterers  when  Pra=  l.E-5 
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Figure  1.33  The  reproduced  figure  for  4-pulsc  slow-lliicliialing  dominant  plus 
Rayleigh  scattercrs  when  Pfj,  =  l.E-8 
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Figure  1.35  The  reproduced  figure  lor  4-pulse  ftisl-fluciUtUing  dominant  plus 
Rayleigh  scatterers  when  Pf.,=  l.E-2 


Figure  1.36  The  reproduced  ligure  for  4-pulse  fast-fluctualing  dominant  plus 
Rayleigh  scatterers  when  Pf.,  =  l.E-5 
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Figure  1.37  The  reproduced  figure  for  4-pulsc  Idsl-fliicuiating  dominant  plus 
Rayleigh  scatterers  when  1  E-S 


Figure  1.38  The  reproduced  figure  for  16-pulse  fasl-lluctualing  dominant  plu 
Rayleigh  scatterers  when  Pm  =  1  .E-2 


Figure  1.39  The  reproduced  figure  for  16-pulse  fasl-llucliiuting  dominant  plus 
Rayleigh  scatterers  when  Pf,,  =  l.E-5 
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Part  II 


THE  EFFECTS  OF  NON-U NIFORMLY  DISTRIBUTED  ENERCJY 


ON  RADAR  DETECTION  PERFORMANCE 


INTRODUCTION 

Radar  detection  peifoimance  has  been  studied  for  several  decades,  e  g.  |11.  In  the 
existing  works,  it  is  generally  assumed  that  the  energy  is  uniformly  distributed  over  the  cells. 
In  1 1 ),  the  probability  of  detection  was  plotted  vs.  signal-  to-noisc  ratio  for  different  number 
of  cells  and  different  target  models  under  the  assumption  of  uniformly-distributed  energy. 
The  target  models  studied  in  ( 1]  includes; 

1)  Single  pulse,  constant  amplitude  scatterers, 

2)  Single  pulse,  Rayleigh  scatterers, 

,1)  Single  pulse,  dominant  plus  Rayleigh  scatterers, 

4)  Multi-pulses,  constant  amplitude  .scatterers, 

.'>)  Multi-pulses,  slow  fluctuating  targets,  Rjiylcigh  scatterers, 

6)  Multi-pulses,  fast  lluctuating  targets,  Rayleigh  scatterers, 

7)  Multi-pulses,  slow  fluctuating  targets,  dominant  plus  Rayleigh  scatterers, 

8)  Multi-pulses,  fast  fluctuating  targets,  dominant  plus  Rayleigh  scatterers. 

However,  in  practical  applications,  the  cells  may  have  different  energy;  i.e.,  the  energy 
is  not  uniformly  distributed  over  the  cells.  In  these  ca,scs,  the  actual  probability  of  detection 
will  differ  from  that  with  uniformly  distributed  energy  and  most  of  the  formula  involved  in 
computing  the  probability  must  be  modified. 

A  difficulty  in  studying  the  effects  of  energy  distribution  is  caused  by  too  many  possible 
distribution  manors  of  the  energy.  To  overcome  this  difficulty,  we  have  to  develop  a  simple 
energy  distribution  model  which  can  well  represent  the  practical  cases. 

In  this  part  of  the  report,  we  will  make  all  the  ncces.sary  modifications  to  the  formula 
involved,  develop  a  simple  and  realistic  model  of  energy  distribution  which  is  described  by  only 
twoparameters.  Then  we  will  compute  the  probability  of  detection  for  different  target  models, 
and  discuss  the  change  of  detection  probability  caused  by  the  energy  distribution.  Finally,  the 
change  will  be  interpreted  by  the  results  corresponding  to  uniformly  distributed  energy. 


FORMULATION 

To  develop  Ihe  energy  disli  ibiilion  model,  we  assume  lhaf  Ihe  total  number  of  resolution 
cells,  M,  is  given  and  the  energy  of  the  cells  is  fallen  into  a  preset  number  of  energy  levels.  The 
energy  of  the  second  level  is  twice  the  energy  of  the  1st  level,  the  energy  of  the  .3rd  level  is  three 
times  the  energy  of  the  1st  level,  and  so  on. 

Generally,  the  energy  of  more  cells  is  in  the  middle  levels,  and  energy  of  only  a  few  cells 
is  in  the  upper  levels  and  lower  levels.  Therefore,  we  assume  that  Ihe  number  of  the  cells  with 
specific  energy  is  normally  distributed  over  energy  levels.  The  number  of  cells  at  energy  level 
El  is  given  by 

=  S  exp  ^ ^ j  (/ =  1,  2, ...,  N^)  (2.1) 

where  N/  is  the  total  number  of  energy  levels.  Since  M/  is  generally  not  an  integer  at  an  energy 

level,  it  has  to  be  rounded.  In  Eq.  (2.1),  E  and  a  arc  free  parameters  to  control  the  shape  of 
the  distribution,  and  S  is  so  selected  that  the  total  number  of  cells  equals  to  the  given  number, 
M.  Figure  2.0  shows  an  example  with  eight  energy  levels. 

Note  the  difference  between  energy  distribution  over  cells  and  the  distribution  of  the 
number  of  cells  over  the  energy  levels.  When  a  is  small,  the  number  of  cells  is  distributed 

nearby  E;  i.e.,  Ihe  number  of  cells  is  unevenly  distributed  over  the  energy  levels.  In  this  case, 

all  the  cells  have  almost  the  same  energy,  E;  i.e.,  the  energy'  is  almost  uniformly  distributed 
over  the  cells. 

When  the  energy  is  not  uniformly  distributed  over  the  cells,  most  of  the  formula  given 
in  Ref.  f  1)  must  be  modified.  The  probability  of  detection  is  given  by 


N, 


Pd=  i-.n(i-p,„,.) 


M 


i=  1 


(2.2) 


where  Mj  is  the  number  of  cells  at  the  ith  energy  level.  Corresponding  to  different  target 


models,  Ihe  formula  (o  compute  Piim.  are  given  by 

(1)  Single  Pulse,  Constant  Amplitude  Scatterers 

^  ''"’i  "  '  -  ’  ^2.3) 

where  Ty(i,  j,  kj  is  the  incomplete  Toronto  function;  No  is  the  total  noise  energy;  E.  is 

Ihe  total  energy  at  ilh  level;  is  Ihe  threshold.  The  formula  to  compute  E.  and  Yh  will  be 
given  later. 

(2)  Single  Pulse,  Rayleigh  Scatterers 


Yu 


to 


(5)  Multi-Pulses,  Slow  Fluctuating  targets,  Rayleigh  Scatterers 


Pdm  =  1  -  I  f  f  l  - 1 

i  ’  )  V  NE/lMjNoJ 

Yb 


N-1 


exp 


(  I  +NEy/ 2M  iN(,  )  *  *  (>/N^I 


(1+2MjN,/NE^ 


-  ,N-2| 

E)  ) 


where  I(x,  j)  is  the  incomplete  gamma  function. 


(6)  Multi-Pulses,  Fast  Fluctuating  targets,  Rayleigh  Scatterers 

Yb 


dm 


(l  +  Ei/2MiN^ 


,N-1  j 


(2.8) 


(7)  Multi-Pulses,  Slow  Fluctuating  targets,  Dominant  plus  Rayleigh  Scatterers 


where 


Pdm.=  f  P(Y)dY 

'  Y„ 

Yfl+  NE./4M.N  ) 
=  V  i  »  0  / 


(2.9a) 


N-2 


P(Y)  =  - 5^ _ ! _ -  .  K  ,  exp^ - ^ 

(l  +  NE./4MjN^^  Ul  +  NE/^MjN^ 

f  __J _  ^N-1 

(N-2)[1-HnEV4M.nJ  exp^ - 


and 


) 


+ 


(l  +  NE./4MjN^ 

N-l  , 

Y  cxp(-Y) 
(N-2)!  (1  +  NE./4MjN^ 


(l  +  NE./4MjN^ 


) 


(2.9b) 


K  =  I 


(l+4MiN^NEj) 


(2.9c) 


(2.7) 
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(8)  Mulli-Pulscs,  Fast  Fluctuating  targets.  Dominant  plus  Rayleigh  Scatterers. 


^  dm .  —  I 


N! 


N 

I 


)!  (  4MiN„  ) 


(l  +  E./4MjN^^  k  =  ()  klCN-k) 

•  I  (  , _ ^ -  ,  N  +  k-1  )  (2 

'^'^NTkd+E.MMjN^)  J 

The  threshold,  Y|„  is  related  to  false  alarm  probability,  Pfn,  by 


10) 


Pf.  = 


for  single  pulse  cases 


(2.11a) 


and 


P,  =  1  - 

fa 


^1-exp  (-Y,)  j  ^ 

,  N-l^  ^  for  multi-pulse  cases  (2  I  lb) 


7'he  total  signal  energy  in  a  given  energy  level,  E„  is  related  to  the  total  signal  energy,  E ,  by 

(2.12) 


-  N,  _  N, 

E=l  E.  =  X'MjEi, 


i=l  '  (=1 

where  Mj  is  the  number  of  cells  in  the  ith  energy  level,  Ei,  is  the  base  energy.  Then 

N, 


E,  =  i  M  jE,,  =  i  M  iE 


7 


i  -  /  2  j  M  j 

i  =  l 


(2.13) 


RluSULTS 

Corresponding  to  different  target  models,  different  number  of  cells,  different  pulse 
numbers,  and  different  energy  levels,  we  have  plotted  a  number  of  figures  which  show  the 
probability  of  detection  vs.  the  ratio  of  total  energy  to  noise.  We  plotted  eight  sets  of  figures 
corresponding  to  eight  target  models.  Four  pulses  arc  assumed  for  mult  i-pulsc  cases.  Each  set 
consists  of  seven  figures,  which  arc  described  as  follows: 
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(a)  Figures  2.1a-2.8a  show  (he  comparison  of  deieclion  performance  helween  uniformly 
dislribuled  energy  and  non-uniformly  disiribuled  energy  when  M  =  2  and  a  false  alarm 
probability  of  0.01.  This  false  alarm  probability  corresponds  to  a  Yb  value  of  .*>.296  for 
single  pulse  cases  and  a  Yh  value  of  10.97  for  miilli-piilse  cases.  M  is  the  total  number  of 
cells  and  Yb  is  the  threshold.  When  energy  is  uniformly  distributed  over  the  cells,  both  cell 
have  the  same  energy;  when  energy  is  non-uniformly  distributed  over  the  cells,  ( he  energy 
of  one  of  the  cell  is  twice  that  of  (he  other  cell. 

(b)  Figures  2.1b-2.8b  show  the  same  ca.ses  as  (hose  in  (a)  except  a  false  alarm  probability  of 

1  .E-8.  This  false  alarm  probability  corresponds  to  a  Yb  value  of  19. 1 1  for  single  pulse  cases 
and  a  Yb  value  of  27.35  for  multi-pulse  ca.ses. 

(c)  Figures  2.1c-2.8cshow  (he  comparison  of  detect  ion  performance  between  different  values 

of  CT  in  Eq.  (2. 1)  when  N/  =  8,  M  =  10,  E=  5,  and  a  false  alarm  probability  of  0.01.  This 
false  alarm  probability  corresponds  to  a  Yb  value  of  6.903  for  single  pulse  cases  and  a  Yb 

value  of  13.05  for  multi-pulse  cases.  E  is  the  energy  level  with  largest  number  of  cells  (see 
Eq.  (2.1)),  N/  is  the  number  of  energy  levels,  M  is  the  total  number  of  cells.  The  energy 
is  more  uniformly  distributed  over  the  cells  when  o  is  smaller. 

(d)  Figures  2.  Id-2.8d  show  the  same  ca.ses  as  those  in  (c)  except  a  false  alarm  probability  of 

l.E-8.  This  false  alarm  probability  corresponds  to  a  Ybvalucof  20.72  for  single  pulse  cases 
and  a  Yi,  value  of  29.15  for  multi-pulse  cases. 

(e)  Figures  2.  le-2.8e  show  (he  comparison  of  detect  ion  performance  between  different  values 

of  Ewhen  N/  =  8,  M=  10,  a  =  0.5  and  a  false  alarm  probability  of  0.01.  This  false  alarm 
probability  corresponds  to  a  Yb  value  of  6.903  for  single  pulse  cases  and  a  Yb  value  of  1 3.05 
for  multi-pulse  cases. 

(I)  Figures  2.1f-2.8f  show  the  same  cases  as  those  in  (c)  except  a=  3.0. 
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(g)  Figures  2.1g-2.8g  show  I  he  same  cases  as  I  hose  in  (e)  except  a  =  2.0  and  a  false  alarm 
probability  of  l.E-8.  This  false  alarm  probability  corresponds  to  a  Y|,  value  of  20.72  for 
single  pulse  cases  and  a  Yi,  value  of  29.15  for  mulli-pul.se  ca.ses. 

In  the  figures,  the  horizontal  axis  denotes  the  ratio  of  total  signal  energy  to  total  noise 

energy  (E/No)  in  decibel.  The  vertical  axis  denotes  the  probability  of  detection  in  percentage. 
Y|,  is  the  threshold,  M  is  the  total  number  of  cells,  “Level”  is  the  total  number  of  energy  levels, 

“Mean”  denotes  E. 

The  following  eight  sets  of  figures  describe  different  target  models; 

Figures  2.1a-2.1g:  single-pulse  constant  amplitude  scatterers; 

Figures  2.1a-2.2g:  single-pulse  Rayleigh  scatterers; 

Figures  2.3a-2..5g;  single-pulse  dominant  plus  Rayleigh  scatterers; 

Figures  2.4a-2.4g;  four-pulses  constant  amplitude  scatterers; 

Figures  2.5a-2.5g:  four-pulse  slow-fluctuating-target  Rayleigh  scatterers; 

Figures  2.6a-2.6g:  four-pulse  fast-fluctuating-targct  Rayleigh  .scatterers; 

Figures  2.7a-2.7g:  four-pulse  slow-fluctuating-target  dominant  plus  Rayleigh 
•scatterers;  Figs.  1.1c  and  2.7f  are  omitted  to  save  computational  efforts. 

Figures  2.8a-2.8g:  four-pulse  fast-flucliiating-targct  dominant  plus  Rayleigh 
scatterers. 

Figs.  2.9-2. 12  show  the  energy  distribution  of  the  cells  over  energy  levels  for  various  o 

and  E  when  the  total  number  of  cells  is  10  and  the  total  number  of  energy  levels  is  8.  When 

(T  is  small,  the  cells  are  distributed  over  a  few  energy  levels  nearby  E,  i.e.,  the  cells  have  nearly 
equal  energy.  This  situation  is  similar  as  uniformly-distributed  energy  over  the  cells.  When 

0  =  0. 1,  all  the  cells  are  in  the  energy  level  E;  i.e.,  the  energy  is  uniformly  distributed  over  the 


Thble  2. 1  shows  Uie  distribution  of  the  cells  over  energy  levels  for  different  values  of  a 


when  total  number  of  cells  is  10  and  total  number  of  energy  level  is  8  and  E  =  5. 


Thble  2.1  Distribution  of  the  10  Cells  over  8  Energy  Levels  for  Different 
Values  of  a  when  E  =  5. 
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DISCUSvSION 

The  figures  in  this  part  of  (he  report  show  the  following  common  phenomenon  for  every 
target  models: 

When  Yj,  is  small,  the  probability  of  detection  is  not  significantly  affected  by  energy 
distribution.  When  Yi,  becomes  larger,  the  effects  of  cncrg>'  distribution  become  more 
significant.  For  example,  when  Y  =  29.15  and  the  total  number  of  cells  is  10,  in  the  ca.se  of 

4-pul.se  constant  amplitude  scattcrers  the  detection  probability  is  about  20%  when  E/No=  15 
db  and  the  energy  is  uniformly  distributed  over  the  cells  (sec  Fig.  2.4d,  (t= 0. 1).  However,  when 
(T=  3.0  (energy  is  not  uniformly  distributed  over  the  cells),  the  detection  probability  becomes 
about  55%. 

In  addition  to  this  common  phenomenon  for  every  target  mode.  We  can  also  find  that 
the  effects  of  energy  distribution  are  also  highly  affected  by  target  models.  We  found  that 
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mulli-  pulse  consfanl  amplitude  scatterers  are  most  seriously  affeeted  by  energy  distribution; 
single  pulse  constant  amplitude  saUterers  are  the  second;  then  it  comes  in  the  order  of 
multi-pulse,  fast-fluctuating-target  dominant  plus  Rayleigh  scatterers,  multi-pulse, 
fast-fluctuating-target  Rayleigh  scalterci^,  multi-pulse  slow-fluctuating-target  dominant 
plus  Rayleigh  scatterers,  multi-pulse  slow-fluctuating-target  Rayleigh  scatterers, 
single-pulse  dominant  plus  Rayleigh  scatterers,  and  finally  single-pulse  Rayleigh  scatterers. 
For  constant  amplitude  scatterers  (single-pulse  or  multi-pulse),  the  effects  of  the  energy 
distribution  are  similar  no  matter  the  probability  of  detection  is  low  or  high.  However,  for 
other  target  models  (Rayleigh  scatterers  and  dominant  plus  Rayleigh  scatterers),  the  effects 
of  energy  distribution  is  significant  in  some  of  the  probability  range  and  less  significant  in  the 
other  ranges.  Generally,  non-uniformly  distributed  energy  will  increase  the  probability  of 
detection  and  the  effects  are  more  significant  in  low  probability  ranges.  The  reason  is  that  the 
threshold  is  hard  to  across  in  that  range  and  the  non-uniformly  distributed  energy  increase  the 
chance  to  across  the  threshold  significantly  for  the  cells  with  higher  energy. 

Let’s  define  the  "cell  probability  of  detection”  as  Pdmj  when  the  energy  is  uniformly 
distributed.  Then  the  effects  of  the  energy  distribution  can  also  be  explained  by  the  shape  of 
the  cell  probability  of  detection.  First,  it  is  helpful  to  consider  the  problem  with  only  two  cells. 
The  probability  of  detect  ion  of  the  first  cell  ispi,und  the  probability  of  detection  of  the  second 
cell  is  p2.  When  the  energy  is  uniformly  distributed,  pi  =  P2  =  P  (cell  probability  of  detection). 
The  total  probability  of  detection  is  given  by: 

P  =  1-(1-P,)(I-P2)  =  2p-p^  (2.14) 

Now,  let’s  assume  the  energy  is  not  uniformly  distributed  and  the  first  cell  has  lower  energy 
and  the  second  cell  has  higher  energy.  Then  pi  becomes  (p-5pi),  and  p2  becomes  (p  +  5p2)- 
The  total  probability  of  detection  becomes 

P  +  5P  =  I-(I-P,)(1-P2)  =  P  +(l-p)(8p2-8p, )  (2.15) 

therefore 
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(2.16) 


8P  =  (l-p)(8p2-8p,) 

We  can  see  the  (rends  of  8P  from  (he  plot  of  (he  p  vs.  signal-(o-noi.  ’  :  lio  in  linear  scales. 
When  (he  cur\'e  is  concavc-up,  8p2  >  8pi,  and  (he  non-  uniformly  distributed  energy  will 
increase  (he  probability  of  detection.  When  the  curve  is  concave-down,  8pi  >  8p2 ,  and  (he 
non-uniformly  distributed  energy  will  decrea.se  the  probability  of  detection.  The  increment 
or  decrement  will  clearly  depend  on  how  large  is  the  curvatu  re  of  t  he  curve  of  p,  how  t  he  energy 
is  distributed  and  how  large  p  is. 

As  an  example,  Figure  2.1.^  plots  the  cell  probability  of  detection  with  different  values 
of  Yb  (6.90.3, 20.72)  for  single-pulse  Rayleigh  scattcrers  when  M  =  10.  The  figure  shows  that 
below  23  db  of  signal-to-noise  ratio  the  curve  is  rclat  ivcly  straight  when  Yb  =  6.903;  therefore, 
the  probability  of  detection  is  not  much  affected  by  the  energy  distribution  (see  Fig.  2.2c).  For 
Yb  =  20.72,  Figure  2. 13  shows  (hat  the  curvature  of  the  curve  is  relatively  large  in  certain  ranges 
of  signal-to-noise  ratio.  The  curve  changes  from  concave-up  to  concave-down  at  about  24 
db.  Therefore,  the  probability  of  detection  is  moresignificantly  affected  by  energy  distribution 
in  low  probability  range  and  is  not  much  affected  by  energy  distribution  nearby  24  db.  And 
when  (he  signal-to-noise  ratio  is  below  24  db,  non-uniformly  distributed  energy  incrca.scs 
probability  of  detection  and  when  the  ratio  is  above  24  db,  non-uniformly  distributed  energy 
decreases  probability  of  detection  (sec  Fig.  2.2d). 


SUMMARY 

In  realistic  radar  detection,  the  cells  have  different  energy.  This  fact  makes  modification 
to  (he  ordinaiy  theory  and  formula  nccc.ssary.  In  this  report,  wc  use  a  realistic  model  of  energy 
distribution  to  investigate  the  probability  of  detection  for  various  target  models. 

We  have  found  a  a)mmon  phenomena  for  every  target  model:  the  effect  of  energy 
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distrihulion  will  become  more  significanl  when  Ihe  Ihreshold,  Y|„  is  increased.  However,  the 
effecl  of  energy  disirihiilion  will  also  largely  depend  on  the  target  models.  Constant  amplitude 
scatterers  are  most  seriously  affected  by  energy  distribution,  dominant  plus  Rayleigh 
scatterers  are  Ihe  next,  and  Rrryleigh  scaflerers  arc  least  affected  by  energy  distributions. 
Multi-pulse  cases  arc  more  seriously  affected  by  energy  distribution  than  single-pulse  cases. 
And  fast  fluctuating  targets  are  more  seriously  affected  by  energy  distribution  than  slow 
fluctuating  targets. 

The  effect  of  energy  distribution  can  be  interpreted  by  a  plot  of  so-called  “cell 
probability  of  detection”  vs.  signal-to- noise  ratio  in  linear  scales.  When  the  curve  is 
concavc-up,  the  non-uniformly  distributed  energy  will  increase  the  probability  of  detection. 
When  the  curve  is  concave-down,  Ihe  non-uniformly  distributed  energy  will  decrease  the 
probability  of  detection.  The  increment  and  decrement  will  depend  on  the  cuivature  of  (he 
curve,  Ihe  energy  distribution  and  Ihe  cell  probability  of  detection  at  the  given  energy  -lo- 
noise  ratio.  In  most  cases,  non-uniformly  distributed  energy  will  increase  the  probability  of 
detection. 
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Figure  2.0  An  Example  of  Cell  Distribution  over  Eight  Energy  Levels 
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Figure  2.1a  Comparison  between  Uniformly  Distributed  Enevg>'  and  Non-Uniformly 

Distributed  Energy  for  Single-Pulse  Constant  Amplitude  Scattcrei's  when  M  =  2  and 
Pf3=I.E-2 
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Figure  2.1b  Comparison  between  Uniformly  Distributed  Encrg>'  and  Non-Uniformly 

Distributed  Energy  for  Single-Pulse  Constant  Amplitude  Scatterers  when  M  =  2  and 
Pfa=l.E-8 
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Figure  2.  Id  Comparison  between  Different  Vrilues  of  a  for  Single-Pulse  Constant  Amplitude 
Scatterers  when  N/  =  8,  M  =  10,  E=  5  and  Pf.,  =  I .E-8 
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Figure  2.  Ic  Comparison  between  Different  Values  of  E  forSingIc-PiiIsc  Constant  Amplitude 
Scatterers  when  N/  =  8,  M  =  10.  ct  =  0.5  and  P,,^  =  I.E-2 
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Probability  of  Detection  I%1 


Figure  2.  If  Comparison  beKvccn  Different  Values  of  E  for  Single-Pulse  Constant  Amplitude 
Scatlerers  when  N/  =  <S,  M  =  10,  ct= 3.0  and  Pf.^  =  l.E-2 


Figure  2.  Ig  Comparison  between  Different  Values  of  E  for  Single-Pulse  Constant  Aniplilude 
Scalterers  when  N/  =  .S,  M  =  JO,  o-=2.0and  Pf.T=  J.E-S 
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Figure  2.2a  Comparison  between  Uniformly  Disiribulcd  Energy  and  Non-Uniformly 

Distributed  Energy  for  Single-Pulse  Rayleigh  Scatterers  when  M  =  2  and  Pf^  =  l.E-2 
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Figure  2.2b  Comparison  between  Uniformly  Distributed  Energy  and  Non-Uniformly 

Distributed  Energy  for  Single-Pulse  Rayleigh  Scatterers  when  M  =  2  and  P(v,  =  l.E-S 
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Figure  2.2c  Comparison  between  Different  Values  of  cr  for  Single -Pulse  Rayleigh  Scaltercrs 
when  N/  =  8,  M  =  10,  E=  5  and  Pf.,  =  I.E-2 
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Figure  2.2d  Comparison  betvveen  Different  Values  of  a  for  Single-Pulse  Rayleigh  Seatterers 
when  N/  =  8,  M  =  10,  E=  5  and  Pf,,  =  I.E-8 
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Figure  2.2e  Comparison  between  Different  Values  of  E  for  Single-Pulse  Rayleigh  Scattercrs 
when  N/  =  8,  M  =  10,  CT=  0.5  and  P(a  =  l.E-2 
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Figure  2.2f  Comparison  between  Different  Values  of  E  for  Singlc-Piilsc  Rayleigh  Scattercrs 
when  N/  =  8,  M=  10,  ct=3.0  and  Pfa=  l.E-2 
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Figure  2.2g  Comparison  between  Different  Values  of  E  for  Single-Pulse  Rayleigh  Scatterers 
when  N/  =  8,  M  =  10,  CT=  2.0  and  Pf^  =  l.E-8 
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Figure  2.3a  Comparison  between  Uniformly  Distributed  Energy  and  Non-Uniformly 
Distributed  Energy  for  Single-Pulse  Dominant  plus  Rayleigh  Scatterers 
when  M  =  2  and  Pf^  =  l.E-2 
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Figure  2.3b  Comparison  between  Uniformly  Distributed  Energy  and  Non-Uniformly 
Distributed  Energy  for  Single-Pulse  Dominant  plus  Rayleigh  Scatterers 
when  M  =  2  and  Pfa  =  I  .E-8 
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Figure  2.3c  Comparison  between  Different  Values  of  u  for  Single-Pulse  Dominant  plus 
Rayleigh  Scatterers  when  N/  =  8,  M=  10,  E=5  and  Pf,T=  I.E-2 
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Figure  2.3d  Comparison  between  Different  Values  of  cr  for  Single-Pulse  Dominant  plus 
Rayleigh  Scatterers  when  N/  =  8,  M  =  10,  E  =  5  and  Pr.,  =  l.E-8 
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Figure  2.3e  Comparison  between  Different  Values  of  E  for  Single-Pulse  Dominant  plus 
Rayleigh  Scatterers  when  N  /  =  8,  M  =  10, 0=  0.5  and  Pfa  =  l.E-2 
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Figure  2.3f  Compari.soii  between  Different  Values  of  E  for  Single-Pulse  Dominant  plus 
Rayleigh  Seatterers  when  N/  =  8,  M  =  It),  (t  =  3.0  and  Pf.,  =  l.E-2 


Figure  2.3g  Comparison  bclwccn  Different  Values  of  E  for  Single-Pulse  Dominant  pi 
Rayleigh  Scatlerers  when  N/  =  8,  M  =  10,  0=  2.0  and  Pfn  =  l.E-8 
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Figure  2.4c  Comparison  bciwccn  Different  Values  of  cr  for  4-Pulsc  Constant  Amplitude 
Scatterers  when  N/  =  8,  M  =  10,  E=5  and  Pf^  =  I.E-2 


Figure  2.4d  Compari-son  between  Different  Values  of  a  for  4-Piilse  Constant  Amplitude 
Scatterers  when  N/  =  8,  M  =  10,  E  =  5  and  Pf;,  =  I.E-8 
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Figure  2.4e  Comparison  between  Different  Values  of  E  for  4-Pulse  Constant  Amplitude 
ScattererswhenN/  =  8,  M  =  10,  (t=0.5  and  ?f„  =  I.E-2 
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Figure  2.4f  Comparison  between  Different  Values  of  E  for  4  Pulse  Constant  Amplitude 
Scatterers  when  N/  =  8,  M  =  10,  ct  =  34)  and  P,,,  =  I.E-2 
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Figure  2.4g  Comparison  belween  DifferenI  Values  of  E  for  4-Piilsc  Constant  Ampliiiide 
Scatterers  when  N/  =  8,  M  =  10,  ct= 2.0  and  Pf.,  =  l.E-8 
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Figus-e  2.5a  Comparison  belween  Uniformly  DislribuJcd  Energy  and  Non-Uniformly 

Distributed  Energy  for  4-Pulse  Slow-Fluctuating  Rayleigh  Scatterers  when  M 
andPfa=l.E-2 
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Figure  2.5b  Comparison  between  Uniformly  Distributed  Energy  and  Non-Uniformly 
Distributed  Energy  for4-Pul.se  Slow-Fluctuating  Rayleigh  Scatterers  when  M  = 
and  Pfa=  l.E-8 
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Figure  2.5d  Comparison  between  Different  Values  of  a  for  4-Piilse  Slow-Fluctuating 
Rayleigh  Scaltercrs  when  N  /  =  8,  M  =  10,  E  =  5  and  Pf^  =  l.E-8 


Figure  2.5e  Comparison  between  Different  Values  of  E  for  4-PuIsc  Slow-Fluctuating 
Rayleigh  Scatterers  when  N/  =  8,  M  =  10,  ct= 0.5  and  Pf„  =  1  .E-2 
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Figure  2.5g  Comparison  between  Different  Values  of  E  for  4-Pulse  Slow-Fluctuating 
Rayleigh  Scatterers  when  N/  =  8,  M  =  10,  ct=  2.0  and  Pf.^  =  l.E-8 
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Figure  2.6a  Comparison  between  Uniformly  Distributed  Energy'  and  Non-Uniformly 
Distributed  Energy  for  4-Pulse  Fast-Fluctuating  Rayleigh  Scatlerers  when  M  =  2 
and  Pfa  =  l.E-2 
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Figure  2.6b  Comparison  between  Uniformly  Distributed  Energy  and  Non-Uniformly 
Distributed  Energy  for  4-Pul.se  Fast-Fluctuating  Rayleigh  Scatterers  when  M  =  2 
and  Pfa  =  l.E-8 
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Figure  2.6d  Comparison  between  Different  Values  of  o  for  4-Pulse  Fast-Fluctuating 


Rayleigh  Scatterers  when  N  /  =  8,  M  =  10,  E  =  5  and  Pf.,  =  l.E-8 
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Figure  2.6e  Comparison  between  Different  Values  of  E  for  4-Pulse  Fast-Fluctuating 
Rayleigh  Scatterers  when  N/  =  8,  M  =  10, 0=  0.5  and  Pfa  =  l.E-2 
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Figure  2.6g  Comparison  between  Different  Values  of  E  for  4-Pulsc  Fast-Fluctuating 
Rayleigh  Scatterers  when  N/  =  8,  M  =  10,  ct=  2.0  and  Pf„  =  l.E-8 


Probability  of  Detection  1%) 


Energy  to  Noise  (dbl 

Figure  2.7a  Comparison  between  Uniformly  Distributed  Energj'  and  Non-Uniformly 

Distributed  Energy  for  4-Pulse  Slow-Fluctuating  Dominant  plus  Rayleigh  Scatterers 
when  M  =  2  and  Pfg  =  I.E-2 
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Figure  2.7b  Comparison  between  Uniformly  Distributed  Energy  and  Non-Uniformly 

Distributed  Energy  for  4-Pulse  Slow-Fluctuating  Dominant  plus  Rayleigh  Scatterers 
when  M  =  2  and  Pf,  =  I.E-8 
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Figure  2.7c  Comparison  between  Different  Values  of  or  for  4-Pulse  Slow-Fluctuating 
Dominant  plus  Rayleigh  Scatterers  when  N/  =  8,  M  =  It),  E  =  5  and  Pfa  =  l.E-2 
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Figure  2.7d  Comparison  between  Different  Values  of  <j  for  4-Pulse  Slow-Fluctuating 
Dominant  plus  Rayleigh  Scatterers  when  N/  =  8,  M  =  10,  E=  5  and  Pfa  =  i.E-8 
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Figure  2.1%  Comparison  between  Different  Values  of  E  for  4-PuIsc  Slow-Fluctuating 
Dominant  plus  Rayleigh  Scatterers  when  N/  =  8,  M  =  10,  a=  2.0  and  Pfa  =  l.E-8 
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Figure  2.8a  Comparison  between  Uniformly  Distributed  Energy  and  Non-Uniformly 

Distributed  Energy  for  4-Pulse  East-Fluctuating  Dominant  plus  Rayleigh  Scatterers 
when  M  =  2  and  Pfg  =  l.E-2 
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Figure  2.8b  Comparison  between  Uniformly  Distributed  Energy  and  Non-Uniformly 

Distributed  Energy  for  4-Pulse  Fast-Fluctuating  Dominant  plus  Rayleigh  Scatterers 
when  M  =  2  and  Pf^  =  l.E-8 
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Figure  2.8c  Comparison  between  Different  Values  of  cr  for  4-Piilse  Fast-Fluctuating 
Dominant  plus  Rayleigh  Scatterers  when  N/  =  8,  M  =  10,  E= 5  and  Pfa  =  l.E-2 


Figure  2.8d  Comparison  between  Different  Values  of  cr  for  4-Pulsc  Fast-Fluctuating 
Dominant  plus  Rayleigh  Scatterers  when  N/  =  8,  M=  10.  E  =  5  and  Pf,  =  I.E-8 
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Figure  2.8e  Comparison  between  Different  Values  of  E  for  4-Pulsc  Fast-Fluctuating 
Dominant  plus  Rayleigh  Scatterers  when  N/  =  8,  M  =  10,  a=  0.5  and  Pfa  =  l.E-2 


Figure  2.8f  Comparison  between  Different  Values  of  E  for  4-Pulse  Fast-Fluctuating 
Dominant  plus  Rayleigh  Scatterers  when  N/  =  8,  M  =  10,  a=  3.0  and  Pfa  =  l.E-2 
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Figure  2.9  Distribution  of  10  Cells  over  8  Energy  Levels  for  Different  Values  of  E 
when  0=0.5 
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Figure  2.^0  Distribution  of  10  Cells  over  8  Energy  Levels  for  Different  Values 
of  E  when  a=L0 
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Figure  2.11  Distribution  of  10  Cells  over  8  Energy  Levels  for  Different  Values 
of  E  when  cr=2.0 
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Figure  2.12  Distribution  of  10  Cells  over  8  Energy  Levels  for  Different  Values 
of  E  when  op=3.0 
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Figure  2.13  Cell  Probability  of  Detection  in  Linear  Scales  for  Single-Pulse 
Rayleigh  Scatterers 
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